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ABSTRACT 

We present Green Bank Telescope (GBT) observations of the 3i2-3i3 (29 GHz) and 4i 3 -4i4 (48 
GHz) transitions of the H2CO molecule toward a sample of 23 well-studied star- forming regions. 
Analysis of the relative intensities of these transitions can be used to reliably measure the densities 
of molecular cores. Adopting kinetic temperatures from the literature, we have employed a Large 
Velocity Gradient (LVG) model to derive the average hydrogen number density [n(H2)] within a 16" 
beam toward each source. Densities in the range of 10 5 5 -10 cm" 3 and ortho-formaldehyde column 
densities per unit line width between 10 135 and 10 14 5 cm" 2 (km s" 1 )" 1 are found for most objects, in 
general agreement with existing measurements. A detailed analysis of the advantages and limitations 
to this densitometry technique is also presented. We find that H 2 CO 3i2-3i3/4i 3 -4i4 densitometry 
proves to be best suited to objects with Tk ^ 100 K, above which the H2CO LVG models become 
relatively independent of kinetic temperature. This study represents the first detection of these H 2 CO 
iVdoublet transitions in all but one object in our sample. The ease with which these transitions were 
detected, coupled with their unique sensitivity to spatial density, make them excellent monitors of 
density in molecular clouds for future experiments. We also report the detection of the 92~8i A" (29 
GHz) transition of CH 3 OH toward 6 sources. 

Subject headings: ISM: clouds — ISM: molecules — stars: formation 



1. INTRODUCTION 

A thorough understanding of star formation is pred- 
icated on a knowledge of the physical conditions that 
surround the process at each of its stages. The determi- 
nation of these properties has proven to be a challenging 
endeavour for the early phases of star formation because 
molecular hydrogen (H2), the primary constituent of de- 
veloping stars, cannot be measured directly. Instead, the 
excitation properties of trace molecules must be used to 
infer the conditions within the greater cloud. Thus, all 
measurements of density in molecular clouds are subject 
to the biases of the chosen tracer. Even more problem- 
atic are virial density estimates determined via inferred 
mass and source size due to the uncertainties inherent in 
both quantities. Formaldehyde (H2CO) is uniquely sen- 
sitive to spatial density and an ideal probe of molecular 
cloud cores for a number of reasons that are discussed in 

m 

The I10-I11 (5 GHz) and 2 u -2i 2 (14 GHz) K- 
doublet transitions of H2CO have previously been 
used to meas ure density in a number of studies of 
galactic (e.g. iHenkel et al.l Il980t iDickel fc Gossl Il987t 
iTurner et all 119891: iZvlka et al.l ll992TT and extragalactic 
(jMangum et all 12008) star formation regions. Our work 
employs an essentially identical strategy using the 3i2- 
3i 3 (29 GHz) and 4 13 -4 14 (48 GHz) transitions, which 
offer a few advantages over their lower excitation coun- 
terparts. Principally, the higher frequencies translate to 
smaller single-dish beam sizes (GBT beam sizes of 26" 
and 16" for the J = 3 and J — 4 transitions, respec- 
tively, as opposed to 153" and 51" for J = 1 and J = 2 
if -doublets). The added spatial resolution is important 
when considering that stars form in spatially compact re- 
gions within molecular clouds and thus large beam sizes 



may serve to dilute the areas of interest. This point is 
compounded by the fact that the lower excitation tran- 
sitions are by their nature often more sensitive to spatial 
densities and kinetic temperatures lower than those of 
interest in studies of star-forming molecular cores. The 
3i2-3i3 and 4i 3 -4i4 transitions are therefore more effi- 
cient probes of spatial density in this context. 

Few studies of the J = 3 and 4 if -doublet transi- 
tions have been made due partially to their relatively low 
intensities and high centimeter-wav e frequencies . The 
J = 3 transition was first detected by Welch (1970) in ab- 
sorption toward the radio continuum source Sgr B2. Sub- 
sequent experiments were conducted toward the bright- 
est source in our samp l e, Orion-KL (Wilson et all ll 9801 : 
iMvers fc Buxton! 119801: iBastien et all 119851) . after which 
exploration of this transition seemingly ends. We could 
find no previous measurements of the 4i3-4i4 transition, 
whose yet higher frequency and weaker intensity pose 
additional observational difficulties. 

Since our ability to detect these transitions and the 
reliability of deriving density measurements from them 
was uncertain, a sample of very bright and well-studied 
objects were chosen (Table [1]). Strong detections of both 
transitions proved to be fairly easy to obtain, requiring 
an average of 17 min of integration time, and the re- 
sulting spatial density measurements (10 55 -10 65 cm" 3 ) 
are consistent with what is known for molecular cores. 
These results are encouraging for the prospect of future 
experiments. 

Details on the utility of H2CO as a high-density probe 
are given in 32 In SHI our observational and calibration 
procedures are presented with a brief discussion of the 
measurement results in $4] Sj5] describes the details of 
our analysis, including a discussion of the Large Velocity 
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TABLE 1 
Source Positions 



Distance V'lsr 



Source 


a(J2000.00) 


5(J2000.00) 


(PC) 


(km s ) 


W3 IRS 4 


02 h 25 m 


30" 


92 


62°06'20'.'70 


1.950 1 


-35 





W3(OH) 


02 


27 


04 


31 


61 52 23.90 


1.950 1 


-50 





L1448 IRS 3B 


03 


25 


36 


33 


30 45 15.00 


232 2 


8 





NGC 1333 IRAS 4A 


03 


29 


10 


52 


31 13 31.00 


235 2 


6 


8 


NGC 1333 IRAS 4B 


03 


29 


12 


00 


31 13 07.80 


235 2 


6 


8 


L1551 IRS 5 


04 


31 


34 


07 


18 08 05.10 


140 3 


6 


4 


Orion-KL 


05 


35 


14 


46 


-05 22 27.50 


418 4 


9 





Orion-N 


05 


35 


14 


45 


-05 21 03.20 


418 4 


9 





Orion-S 


uo 




13 


55 


-05 24 08.30 


418 4 


6 


5 


OMC-2 IRS 4 


05 


35 


27 


32 


-05 09 39.6 


430 5 


11 





IRAS 05338-0624 


05 


36 


18 


42 


-06 22 06.20 


480 6 


7 


1 


NGC 2024 


05 


41 


44 


53 


-01 55 02.90 


415 s 


11 





NGC 2024 FIR 5 


05 


41 


44 


17 


-01 55 42.50 


415 5 


11 





NGC 2024 FIR 6 


05 


11 


45 


16 


-01 56 04.80 


415 s 


11 





NGC 2071 IR 


05 


17 


04 


85 


00 21 43.00 


390 5 


10 





S255N 


0(5 


12 


53 


67 


18 00 26.60 


1,590 7 


8 





G34.26+0.15 


18 


53 


18 


60 


01 14 58.40 


3,700 s 


58 


7 


S68N 


18 


29 


47 


91 


01 16 46.50 


230 9 


8 


8 


W51M 


19 


23 


11 


01 


14 30 34.10 


5,410 10 


55 





DR 21(OH) 


20 


39 


01 


09 


42 22 48.90 


2,000 n 


-5 





Cep A East 


22 


56 


18 


15 


62 01 46.00 


700 1 


-10 


1 


NGC 7538 IRS 1 


23 


13 


45 


34 


61 28 10.50 


2.650 1 


-57 


3 


NGC 7538 IRS 9 


23 


11 


01 


66 


61 27 20.00 


2.650 1 


-57 


3 



Re fere nces. — ( HReid et al.1 120091 ^IHirota et al. | |2008| f3lKenvon et aT 
IT9911 (4'IKim et al.l 1 2003 (5' Anthon y-Twarod 119821 fejchen et al.l 11995 
(7;Rve1 ct al. 2010 (S^Kuchar & Bania 1994 (O^Eiroa ct al. 2008 (10;Sato ot al 
120101 fll lOdenwald fc Schwartidll993l 



Gradient (LVG) and Local Thermodynamic Equilibrium 
(LTE) approximations used, as well as source-by-source 
comparisons to past studies. A detailed discussion of the 
limitations of measuring density with the J = 3 and 4 
K -doublet transitions of H2CO is provided in EJSJ 

2. FORMALDEHYDE AS A HIGH-DENSITY PROBE 

Formaldehyde (H2CO) is a proven tracer of the high- 
density environs of molecular clouds. It is ubiquitous; 
H2 CO is associated with 80% of the HII regions surveyed 
by iDownes et al.1 (|1980h and possesses a large number 
of observationally accessible transitions from centimeter 
to far-infrared wavelengths. Because H2CO is a slightly 
asymmetric rotor molecule, most rotational energy lev- 
els are split by this asymmetry into two energy levels. 
Therefore, the energy levels must be designated by a to- 
tal angular momentum quantum number, J, the projec- 
tion of J along the symmetry axis for a limiting prolate 
symmetric top, K_i, and the projection of J along the 
symmetry axis for a limiting oblate symmetric top, K + \. 
This splitting leads to two basic types of transitions: the 
high-frequency A J = 1, AiC_i = 0, AK +1 = -1 "P- 
branch transitions" and the lower frequency A J = 0, 
A_ftT_i = 0, AK+i = ±1 "Q-branch" transitions, popu- 
larly k nown as the "if-doublet" t ransitions (see discus- 
sion in lMangum fc Woottenlfl993| ). The P-branch tran- 
sitions are only seen in emission in regions where n(fi2) 
> 10 4 cm -3 . The excitation of the i-T-doublet transitions, 
however, is not so simple. For n(fi2) < 10 5,5 cm" 3 , the 
lower energy states of the I10-I11 through 5i4-5i5 K- 
doublet transitions be come overpopulated due to a colli- 
sional selection effect ()Evans et al.lll975t Garriso n et al.l 
1975). This overpopulation cools the J < 5 if-doublets 



to excitation temperatures lower than that of the cosmic 
microwave background, causing them to appear in ab- 
sorption. For n(H2) > 10 5 5 cm" 3 , this collisional pump 
is quenched and the J < 5 i^-doublets are then seen in 
emission over a wide range of kinetic temperatures and 
abundances. 

3. OBSERVATIONS 

The measurements reported here were made using the 
National Radio Astronomy Observatory (NRACQ) Green 
Bank Telescope (GBT) during the periods 2007 January 
20-31, February 19, and 2008 July 2. Single pointing 
measurements were obtained of the 3i2-3i3(28. 974805 
GHz, 6 B = 26") and 4 13 -4 14 (48.284547 GHz, 6 B = 16") 
if-doublet transitions of H2CO toward a sample of 23 
galactic star-forming regions (Table [T]) . The 29 and 48 
GHz measurements were made using the dual-beam Ka- 
(beam separation 78") and Q-band (beam separation 
58") receivers, respectively, over 50 MHz of bandwidth 
sampled by 16,384 channels. The position-switching 
technique was employed with reference position located 
30' west in azimuth from each source position. The corre- 
lator configuration produced a spectral channel width of 
3.052 kHz, which is approximately 0.03 and 0.02 km s" 1 
at 29 and 48 GHz, respectively. The 9 2 -8i A" (28.969942 
GHz) vibrational transition of CH3OH was also captured 
in the Ka-band spectra. 

To calibrate the intensity scale of our measurements, 
several corrections need to be considered: 

1 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Opacity-Historical opacity estimates based on atmo- 
spheric model calculations using ambient pressure, tem- 
perature, and relative humidity measurements indicated 
that ro at 29 and 48 GHz were, respectively, ~ 0.05 and 
0.25 during our observations. The respective opacity cor- 
rections exp[roCSc(£'L)] for the range of source elevations 
(28-75°) averaged 1.069 and 1.485, and were applied uni- 
formly since the variance due to elevation was less than 
the absolute uncertainty in amplitude calibration (see 
below) . 

Flux.- Assuming point-source emission, one can use 
the current relation (derived from measurement) for the 
aperture efficiency 

T} A = 0.71exp(-[0.0163i/(GHz)] 2 ) 

to convert antenna temperature to flux density. At 28.97 
and 48.28 GHz this yields r\ A = 0.568 and 0.382, re- 
spectively. For elevation 90° and zero atmospheric opac- 
ity, T A /S = 2.846r/ A = 1.61 and 1.09 for 29 and 48 
GHz, respectively. These are the degrees Kelvin per 
Jansky calibration factors used to convert our spectra 
to flux assuming point-source emission. Incorporating 
atmospheric opacity and telescope efficiency, we have 
T% = T A exp(AT )/rn = 1MT A (29 GHz) and 1.50T A 
(48 GHz), where rji = 0.99 for the GBT. Finally, using 
rj m b ~ 1.32r] A} we can write the main beam brightness 



temperature as T mb ~ T* A jr\ mh = 1.45Ta (29 GHz) and 
2.9QTU (48 GHz). 

Source Structure-Since the beam sizes are about 63% 
different, and we anticipate source sizes comparable to 
our beam sizes, it is necessary to scale both measure- 
ments to the same beam size to properly intercompare 
them. The beam coupling correction factor for a Gaus- 
sian source and a Gaussian beam to a measurement 
which has already been corrected to the main-beam effi- 
ciency scale is: 



^source * p bcam 

couple — ofl2 



Therefore, to scale the J = 3 T m b measurements 
(0beam = 26") to a source size equal to the beam size 
at J = 4 (#bcam — 16"), an additional factor of 1.82 
must be applied to the J = 3 T m b measurements. With 
this, we have assumed that each object in our sample 
spans 16". Refer to H6.3I for additional justification of 
this assumption. 

Absolute amplitude calibration.-The GBT absolute 
amplitude calibration is reported to be 10-15% at all fre- 
quencies, limited mainly by temporal drifts in the noise 
diodes used as absolute amplitude calibration standards. 



TABLE 2 
H2CO Measurement Results 









rp* a 
A 


^LSR 


FWHM 


log rAf(ortho-H 2 CO)l b 


Source 


Transition 


(K) 


(km s" 1 ) 


(km s" 1 ) 


[cm -2 (km s" 1 )" 1 ] 


W3 IRS 4 


4l3- 


-4i4 


< 0.090 (3<t) 






,3 17 +0.16 




3l2- 


-3l3 


-0.023 ± 0.010 


-47.28 ± 0.44 


6.66 ± 0.91 








-0.068 ± 0.010 


-42.53 ± 0.05 


1.00 ± 0.12 


10. oo_ 007 


W3(OH) 


4i 3 - 


■4l4 


0.216 ± 0.028 


-47.56 ± 0.15 


6.23 ± 0.15 


i4.is±S:8i 








-0.136 ± 0.028 


-44.46 ± 0.15 


1.46 ± 0.15 


13.98lg;g 




3l2- 


-3l3 


0.124 ± 0.014 


-47.60 c 


6.20 c 


14.09±0.05 








-0.284 ± 0.014 


-44.50 c 


1.50 c 


14.45±0.02 








-0.167 ± 0.014 


-46.96 ± 0.08 


3.20 ± 0.34 


14.22 t° i 


L1448 IRS 3B 


4l3- 


-4i4 


< 0.202 (3<x) 








3l2- 


-3l3 


0.041 ± 0.005 


4.85 ± 0.06 


1.76 ± 0.13 


i3-5i±°;| 


NGC 1333 IRAS 4A 


4l3- 


-4l4 


0.084 ± 0.031 


7.00 c 


2.06 c 


13.481°;" 








0.039 ± 0.031 


7.00 c 


8.50 c 






3l2- 


-3l3 


0.162 ± 0.006 


6.90 c 


2.06 ± 0.06 


13.85±0.02 








0.033 ± 0.006 


6.90 c 


8.50 c 


13 16+ 007 








-0.118 ± 0.006 


7.50 c 


2.00 c 


13.71±0.02 


NGC 1333 IRAS 4B 


4l3- 


-4l4 


0.133 ± 0.029 


7.08 ± 0.09 


2.25 ± 0.26 


13.83+°-°? 




3l2- 


-3l3 


0.036 ± 0.009 


6.90 c 


2.20 c 


13.401^° 


L1551 IRS 5 


4l3" 


-4l4 


0.103 ± 0.039 


6.51 ± 0.07 


0.48 ± 0.14 


13.82+g;^ 




3l2- 


-3l3 


0.050 ± 0.006 


6.34 ± 0.03 


0.73 ± 0.08 


13.65±0.05 


Orion-KL 


4l3- 


-4i4 


1.748 ± 0.064 


10.27 ± 0.02 


1.72 ± 0.04 


15.19±0.02 








0.551 ± 0.064 


6.71 ± 0.14 


13.53 ± 0.46 


15.12±0.05 








0.938 ± 0.064 


8.30 c 


2.71 ± 0.17 


14.92±0.03 




3l2- 


-3l3 


0.670 ± 0.016 


10.41 ± 0.001 


1.62 ± 0.003 


14.93±0.01 








0.563 ± 0.016 


6.72 ± 0.03 


9.58 ± 0.004 


15.32±0.01 








3.972 ± 0.016 


8.16 c 


2.46 ± 0.003 


15.71±0.01 


Orion-N 


3l2- 


-3l3 


0.259 ± 0.012 


9.37 ± 0.02 


1.47 ± 0.05 


13.95±0.02 


Orion-S 


4l3- 


-4l4 


0.793 ± 0.035 


6.57 c 


2.68 ± 0.09 


14.70±0.02 








0.243 ± 0.035 


6.90 c 


6.35 ± 0.40 


14 ig+0.06 




3l2- 


-3l3 


0.575 ± 0.011 


6.39 c 0.00 


2.65 ± 0.04 


14.70±0.01 








0.143 ± 0.011 


6.72 c 


6.35 c 


14 i +0.03 


OMC-2 IRS 4 


4l3- 


-4l4 


0.175 ± 0.026 


11.57 ± 0.05 


1.66 ± 0.13 


13 98+ - 06 




3l2- 


-3l3 


0.060 ± 0.008 


11.50 ± 0.12 


1.41 ± 0.05 


1 3 fi^+O- 05 
1>S.DO_ 06 








0.253 ± 0.008 


11.36 ± 0.01 


6.06 ± 0.49 


14.27±0.01 


IRAS 05338-0624 


4l3- 


-4l4 


0.030 ± 0.029 


9.60 c 


2.48 c 


13 26+ ' 29 
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TABLE 2 — Continued 







n->* a 
A 


^LSR 


FWHM 


log lO""" ri 2 K - jKJ J b 


Source 


Transition 


(K) 


(km s" 1 ) 


(km s ) 


[cm -2 (km s" 1 )" 1 ] 






0.075 ± 0.029 


7.08 c 


1.94 c 


13.65_°;« 




3l2-3l3 


0.019 ± 0.008 


9.60 c 


2.48 c 


13 21+015 






0.070 ± 0.008 


7.08 c 


1.94 c 


13.77±0.05 


NGC 2024 


4l3"4l4 


0.149 ± 0.038 


11.50 ± 0.08 


1.25 ± 0.19 


, _ _ __i_n i n 

13.95_°;i° 


NGC 2024 FIR 5 


4l3"4l4 


0.202 ± 0.041 


11.20 c 


1.97 ± 0.19 


14.34±g;°o 




3l2-3l3 


0.224 ± 0.008 


11.10 c 


2.04 ± 0.05 


14.55±0.01 






-0.062 ± 0.008 


9.10° 


0.58 ± 0.06 


13.99±0.06 


NGC 2024 FIR 6 


4l3-4l4 


0.108 ± 0.030 


11.12 ± 0.12 


2.56 ± 0.27 


13.55_Q.14 




3l2"3l3 


0.208 ± 0.011 


11.10 c 


2.14 ± 0.07 


13.89±0.02 






-0.044 ± 0.011 


9.11 ± 0.12 


0.57 ± 0.11 


1 Q 91 T-U.-LU 

io.z±_ 13 


NGC 2071 IR 


4i3-4i4 


0.097 ± 0.040 


8.33 c 


4.79 ± 1.00 


13.69_g_!| 






0.155 ± 0.040 


9.81 ± 0.13 


1.70 ± 0.39 


13.90_°;i° 




3l2-3l3 


n r\o r I n nno 

0.085 ± 0.008 


o oo I n i o 

8.33 ± 0.18 


A TO I n OT 

4.73 ± 0.27 


1 o tt I n n a 

13.77±0.04 






0.175 ± 0.008 


9.58 ± 0.03 


1.84 ± 0.12 


14 08±0 02 






-0.043 ± 0.008 


16.17 ± 0.09 


1.94 ± 0.26 


13.48 +0 ,'°j; 


DZOOlN 


A „ A 

413-414 


n iQi i n nifl 
u.ioi ± u.uiy 


0.0/ it U.Uo 


o O/i i n 1 o 
Z.y4 it U. lz 


1 Q.9 + - 04 

±o.»z_ Q5 




3l2"3l3 


n 1 no _i_ n nns 
U.lyz ± U.UUo 


q 7o _l n no 
o. 7 y ± U.U2 


o no _i_ n c\a. 
2.x) 2 ± U.Uu 


1 /i n^;_i_n no 
14.U0±U.U2 


Vjo4.zD-hU. 10 


4l3-4l4 


U.14< it U.Uoo 


c;t 70 _i_ no/i 
Oi./i it U.z4 


^ i k _i_ n cq 
D.1D it U.Oo 


14.2U_n.ig 


CROAT 


4l3-4l4 


n no j. n n/in 
u.ny ± u.U4u 


/ .50 it u.iy 


4.1o it U.4Z 


co+0.13 
13.58_ lg 


VV D11V1 


4l3-4l4 


n ot;i -t n non 
U.zDl it U.UoU 


Kfi A A -I- n Ot; 
00.44 it U.ZO 


o.iO it U.OO 


14 4n+ ' 12 

14.4U_ 17 


UK 21((JHJ 


4i3-4i4 


n oon I n no ,1 

0.320 ± 0.034 


a oo I n n*" 

-4.22 ± 0.06 


o no I r\ 1 a 

3.02 ± 0.14 


14 - 54 Io.05 






0.345 ± 0.034 


-1.28 c 


3.40 c 


- . ^ .4-0 04 




3l2"3l3 


0.247 ± 0.005 


-4.45 ± 0.04 


2.89 ± 0.06 


14.59±0.01 






n 1 on _l_ n nnec 
0.182 ± 0.005 


i oo _L n c\a. 
-1.28 ± 0.06 


q .in _l_ n i 1 
3.40 ± 0.11 


1 q vt; J_n m 
lo. (0±U.U1 


Cep A East 


4 13 -4i4 


0.078 ± 0.033 


-10.40 c 


2.40 c 


13 50+° if 

j-o.<ju_Q 24 






0.042 ± 0.033 


-10.28 ± 0.59 


6.00 c 


1 1 oq+0.25 
io.zo_ 66 




3l2"3l3 


0.036 ± 0.007 


-10.43 ± 0.10 


2.43 ± 0.23 


,0 97 +0.07 
io.zi_ Q og 


NGC 7538 IRS 1 


4l3"4l4 


0.111 ± 0.028 


-57.70 ± 0.14 


4.18 ± 0.36 


14 25+ - 10 




3l2"3l3 


-0.102 ± 0.009 


-59.00 c 


2.19 ± 0.13 


14.39±0.04 






0.110 ± 0.009 


-57.65 c 


4.00 c 


14 43+0.03 


NGC 7538 IRS 9 


4l3-4l4 


< 0.133 (3a) 










0.048 ± 0.009 


-56.82 ± 0.12 


3.31 ± 0.29 


13 39+ 07 



Source radiation temperature corrected for telescope efficiency and 
atmospheric attenuation. Main beam efficiency and beam dilution 
calibrations applied for analysis (sec 
k LTE approximation. Sec ^5.2l for details. 
c Quantity fixed in Gaussian fit. 

4. RESULTS 

Both transitions were detected in 15 of the 23 sources 
in our sample. Three objects yielded nondetections in 
the J = 4 transition (W3 IRS 4, L1448 IRS 3B, and 
NGC 7538 IRS 9) and another five sources were ob- 
served in only one of the two transitions due to time 
constraints. Data reduction was accomplished using 
the CLASS (Continuum and Line Analysis Single-dish 
Software) package from GILDA _. The following spec- 
tra have been smoothed to 0.2526 km s" 1 (3i2-3i3) and 
0.1516 km s" 1 (4i3-4i4) to increase the signal-to-noise ra- 
tio of the individual channels. Each line profile was fitted 
with 1-3 Gaussian components, and each velocity com- 
ponent was checked for consistency with previous mea- 
surements (see §5.4|) . 

Table _ lists the peak intensity (T|), velocity of the 
local standard of rest (Vlsr), and velocity width at the 
half maximum (FWHM) determined via Gaussian fitting 

2 http://www.iram.fr/IRAMFR/GILDAS/ 



for each object. LTE column densities are also listed (see 
£|5.2p . Note that the peak intensities quoted in Table _ 
are given in the T\ scale in the event that a reader may 
wish to apply calibration factors different than those used 
in our analysis (see SJH) . The spectra for sources detected 
in both the 3i2-3i3 and 4i3-4i4 transitions are displayed 
in Figure _ while spectra from sources detected or ob- 
served in only one of the two transitions are shown in 
Figure _ 

The 92-8i A" transition of CH3OH was detected in 6 
sources. Since no further analysis is presented, Gaussian 
components were not fit to each spectrum and the basic 
line parameters, determined via direct measurement of 
the line peak and width, are listed in Table _ with the 
spectra displayed in Figure _ Maser emission of this 
transition was detected toward one source, W3(OH). 

5. ANALYSIS 

5.1. Spatial and Column Density Derivation from LVG 

Modelling 

To derive the spatial density [H_ number density, 
rj(Il2)] and ortho-H2CO column density in our sample 
of galactic star-forming regions, we have used a model 
which incorpor ates the Large Velocity Gradient (LVG) 
approximation (|Sobolevl 119601 ) to radiative transfer in 
molecular clouds. The detailed properties of our imple- 
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Fig. 1. — Spectra for sources toward which both transitions were detected. 
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Fig. 2. — Spectra for sources toward which only one transition 
was detected or observed. 
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Fig. 3. — Methanol spectra. 



mentation of the LVG app roximation are described in 
iMangum &; Woottenl ([1993) . In short, the model is used 
to predict the individual brightness temperatures of both 
the J = 3 and 4 transitions as well as the "transition ra- 
tio," Rj = / T mb (3i2 - 3i 3 )di// / T mb (4i 3 - i 14 )di/, for a 
given set of physical conditions: spatial density, column 
density, kinetic temperature, and optical depth. The so- 
lution is then indicated by the intersection of the bright- 
ness temperature and Ri predictions. This is illustrated 
by Figure [U which is discussed in more detail at the end 
of this section. It is important to note the behavior of 
the transition ratio. Since we calculate Ri as (lower exci- 
tation) /(higher excitation), a higher ratio implies lower 
spatial density and /o r kinet ic temperature. As noted by 
IMangum &; Woottenl (fl993). one of the major sources of 
uncertainty in an LVG model prediction of physical con- 
ditions is the uncertainty associated with the collisional 
excitation rates used. 

Our model incorporates the scaled excitat i on rat es in- 
volving H 2 CO and He calculated by iGreer] (|1991|) . An 
important point regarding our implementation of the 
LVG model is the scaling of the calculated H 2 CO/Hc 
excitation rates to those appropriate for coll i sions with 
H2. Following the recommendation of iGreenl (|1991l ). we 
scale the calculated He rates by a factor of 2.2 to account 
for (1) the reduced collision velocity of He relative to H 2 , 
which scales as the inverse-square-root of the masses of 
He and H 2 , and (2) the la rger cross section of the H 2 
molecule (~ 1. 6: lNerjHl975f) relative to He. With these 
scaling factors. fGreenl |l991i) suggests that the total col- 
lisional excitation rate for a given H 2 CO transition is 
accurate to ~ 20%. Thus the physical conditions pre- 
dicted by our LVG model are limited to an accuracy of 
no better than 20%. It should be noted that the col- 
lision rates of H 2 CO with H 2 have been rederiv e d wit h 
a claimed accuracy of 10% by iTroscompt eTall ([200l . 
though only for temperatures < 30 K. 

It is also important to note that the LVG model for 
the 3i 2 -3i3 and 4i3-4i4 transitions of H 2 CO is somewhat 
dependent on kinetic temperature (Tk). As such, it was 
important in this study to carefully select temperature 
estimates from the literature that best reflect the mate- 
rial being traced by our observations. The temperatures 
used can be found with their object references in H5A{ 
The Tk dependence of the LVG model and its effect on 
our measurements is explored in detail in §6.11 

LVG results for the 15 sources toward which both tran- 
sitions were detected are displayed in Table |4] A range 
of derived densities based on the range in our adopted 
kinetic temperatures is provided in column 5 to illustrate 
the kinetic temperature dependence of the LVG model. 
In absence of error bars on the Tk values assumed from 
the literature, a range of ± 50% was used. Column 6 
includes the range of density estimates found in selected 
literature (see t j5.4p . Note that the uncertainties asso- 
ciated with our results reflect only measurement error. 
Additional uncertainties associated with our kinetic tem- 
perature and spatial extent assumptions are described in 

m 

For 6 of 21 velocity components, the observed transi- 
tion ratio [Rj = J T mh (3i 2 - 3i 3 )g?j// J T mfc (4i 3 - 4u)di/] 
and brightness temperatures (T m b) were not consistent 
with LVG predictions and a Local Thermodynamic Equi- 
librium (LTE) approximation was employed to place a 
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limit on the density (see ij5.2[) . In general, this occurred 
when the was low (< 1), meaning that the higher ex- 
citation J — A transition was unexpectedly observed to 
be significantly brighter than the J — 3. Figure 0] dis- 
plays |LVG Prediction — Observation! for the transition 
ratio and J = 4 brightness temperature in units of the 
measurement error a for two sources. OMC-2 IRS 4 (left 
panel) is indicative of those sources with observed Ri ~ 



1-2 for which the density could be fully constrained. The 
solution is indicated by the intersection of the violet and 
shaded (< la) regions. Cep A East (right panel) is repre- 
sentative of those objects with ratios < 1. Note how the 
color and solid-line contours are parallel. Further discus- 
sion of the reasons for these failed model fits is given in 



CH 3 OH 9 2 -8; 



TABLE 3 

A" Measurement Results 





T* 




Vhsn 


FWZI 




Source 


(K) 




(km 


s" 




(km s" 


- 1 ) 


(K km 


r 1 ) 


W3(OH) a 


19.211 ± 


0.014 


-43.40 


± 


0.01 


3.06 


± 


0.01 


12.093 ± 


0.012 


Orion-KL a 


2.483 ± 


0.016 


7.53 


± 


0.01 


18.18 


± 


0.04 


11.089 ± 


0.034 


Orion-S 


0.025 ± 


0.011 


2.39 


± 


0.40 


12.63 


± 


0.82 


0.104 ± 


0.020 


IRAS 05338-0624 


0.021 ± 


0.008 


7.19 


± 


0.36 


9.73 


± 


0.91 


0.061 ± 


0.013 


DR21(OH) 


0.036 ± 


0.005 


-3.28 


± 


0.09 


9.22 


± 


0.23 


0.115 ± 


0.008 


NGC 7538 IRS 9 a 


0.292 ± 


0.009 


-56.20 


± 


0.02 


8.82 


± 


0.06 


0.977 ± 


0.013 



a Multiple velocity components present. See Figurc[3]for spectra. 

5.2. LTE Approximation 

To supplement the LVG model and provide a use- 
ful check on its results, a Local Thermodynamic Equi- 
librium (LTE) approximation to the column density, 
7V(ortho — H2CO)/Ai^ cm" 2 (km s has been calcu- 
lated for each velocity component detected (listed in Ta- 
ble[U): 



QuQ 'K Qnuclear \ ^?ex / 

<3 ™'=K^)" 2r3/2 -"« 4T 

where the line strength S = (J = 3), ^ (J = 4); 
dipole moment /i = 2.331 debye; rotational degeneracy 
9u = 7 (J = 3), 9 (J = 4); K degeneracy g K = 2 (for 
K ^ in symmetric top molecules); nuclear spin degen- 
eracy gnuciear = 3 (for ortho-H^CO); level energy above 
ground E u = 33.479 K ( J = 3), 47.928 K ( J = 4). Note 
that the peak intensity (T m &) is used instead of integrated 
intensity to give column density per unit line width for 
consistency with the LVG model results. Note also that 
the LTE approximation assumes optically thin emission 
with excitation temperature equivalent to kinetic tem- 
perature. 

The former assumption is reasonable given that H2CO 
transitions rarely become optically thick due princi- 
pally to fairly low abundances. This property, in part, 
makes H2CO an ideal probe of the dense regions within 
molecular clouds (see Given the average physical 

conditions derived for our sample [n(H 2 ) ~ 10 6 cm" 3 , 
V(ortho - H 2 CO)/Az/ - 10 14 cm" 2 (km s" 1 )" 1 , T K = 100 
K], the LVG-predicted optical depth (r) ~ 0.078. The 
maximum value found for our sample is 0.27, and tran- 
sitions can generally be considered optically thin for 
t < 0.4. 

The second assumption that T cx = Tk is 
more troublesome. For nQS.2) ~ 10 6 cm" 3 and 



V(ortho-H 2 CO)/A^ ~ 10 14 cm" 2 (km s" 1 )" 1 , the 
excitation temperatures for the J = 3 and 4 JT-doublets 
are on order ^-K, while T ox ~ for the A J = 1, 

AK-i = 0, AK + i = —1 P-branch transitions. There- 
fore, since the ortho-H 2 CO column density is mainly 
dependent upon the excitation temperature of the 
P-branch transitions, T ox < Tk for our sample and 
the resulting LTE column densities must then, strictly 
speaking, be considered upper limits given that T" m b 
oc T cx (l - exp(-r)). Nevertheless, the LTE-predictcd 
column densities are generally within a factor of 2 of 
the LVG results (see ^5 . 31 for more details) . We consider 
the LTE approximation to be a reasonable one in cases 
where the LVG model failed to fully constrain the 
physical conditions. See §5.41 for kinetic temperature 
assumptions (alternatively, Tables HHHJ) ■ 

For those cases (described in the previous section) in 
which LVG modeling was unsuccessful, the unique prop- 
erties of the i^-doublet transitions discussed in §21 were 
used to place a limit on the spatial density. A column 
density was assumed using the LTE approximation de- 
scribed above, and the LVG model was then employed to 
determine the minimum density required for both tran- 
sitions to be observed in emission. 

The LTE approximation was used in a similar manner 
for the 12 velocity components detected in only one of 
the two observed transitions. With the exception of Cep 
A East, these components were detected in the 3i2-3i3 
transition and not the 4i3-4i4. If seen in absorption at 
J = 3, the undetected J = 4 counterpart must also ap- 
pear in absorption because of the excitation of the H2CO 
molecule. If observed in emission at J = 3, the J = 4 
transition may appear in either state because the 4i3-4i4 
transition is subject to absorption beginning at slightly 
higher densities (see M6.2l for further detail). The range of 
densities permitting the J = 3 transition to be in emis- 
sion with the J = 4 in absorption is very small, typically 
spanning 0.2 on the log scale. Rather than dictate this 
range for each object, we have opted to list just the limit 
corresponding to the J = 4 emission/ absorption bound- 
ary. Table [5] displays these results. 



5.3. Singe Transition H2CO Column Density Estimates 

Observations of only one transition were performed to- 
ward 5 sources. For these objects, the LVG model can 
be used to estimate the ortho-H^CO column density by 
assuming a kinetic temperature and spatial density. Ki- 



netic temperatures were culled from the literature and 
can be found with their object references in ^5 .41 A spa- 
tial density of 10 615 cm" 3 was assumed based on an aver- 
age of the results from those sources for which a full LVG 
analysis could be made. These results are presented in 
conjunction with the LTE approximation for comparison 
in Table M 



TABLE 4 

LVG Results: Velocity Components Detected in Both Transitions 















10g[n(H 2 ,)] 






/ T mb (3i2 - 3i3 


)dv 


T K a 




(cm" 3 ) 




Source 


jT mb (4i 3 -4i4 


)du 


( K ) 


^Bcst 


^Rangc 


Literature 13 


W3(OH)= 

\ /emission 




n 73-1-n 03 




110+L 


>4.92+°'jg 

— — U.UO 


>4.88-5.00 


5.95 


W3(OH) a t, sor ptioii 




2.70±0.33 




60±20 


4.11±0.05 


3.97-4.21 


4.70 


NGC 1333 IRAS 4A d 




1.53±0.15 




50 


6.20+°^ 


5.87-7.07 


5.48-6.30 


NGC 1333 IRAS 4B C 




0.34±0.06 




80 


>5.3llg-» 


>5. 10-5.62 


5.48-6.30 


L1551 IRS 5 




0.93±0.33 




10011° 


6.87T>£ 3 


6.49->8.00 


5.78-7.04 


Orion-KL hot core 




0.92±0.02 




300-70 


^ 07+0.10 

°- 9 '-0.09 


6.08 


5.70-6.95 


0ri0n - KL compact ridge 




4.85±0.12 




135±?1 






5.50-5.70 


Orion-KL C , , , . 

extended ridge 




0.45±0.01 




135±?S 


7 01+0.79 
' -0.40 




> 5.00 


Orion-Sp Cak 




0.91±0.02 




100 


6.43->8.00 


6.60-7.78 


Orion-S c . 

wing 




0.74±0.03 




100 


> 4 Q 6 + 02 


>4.66-5.26 




OMC-2 IRS 4 




1.55±0.12 




B0 -10 


5.51±0.09 


5.37-5.59 


5.65-6.28 


IRAS 05338-0624 d 




1.05±0.16 




95±55 


6 4Q+> 151 

°-^ M -0.30 


6.20->8.00 


5.30-6.00 


NGC 2024 FIR 5 




1.45±0.13 




160+ 90 

1DU -60 


°** i -g.l2 


5.31-5.57 


5.95-6.00 


NGC 2024 FIR 6 




2.03±0.23 




40(<160) 


c- oo+O.ll 


4.94 


5.95-6.30 


NGC 2071 IRp eak 




1.54±0.19 




80+ 2 ° 


c 71 +0.18 


5.57-5.85 


5.48-6.48 


NGC 2071 IR w i n g 




1.09±0.13 




80+ 20 . 




6.34-6.73 


> 5.00 


S255N 




1.33±0.06 




70±30 


6 02+ ' 09 


5.78-6.58 


6.48-7.20 


DR 21(OH)mmi 




0.93±0.04 




160 +140 


6 34+O I 6 


6.04->8.00 


5.80-6.10 


DR 21(OH)^ M2 




0.67±0.02 




30(<60) 


>5.58±0.08 


>5.36 


7.00-7.83 


Cep A East c 




0.58±0.11 




' U -30 


>5 42+ 02 

^3-^_0.04 


>5.30-5.91 


6.30-6.78 


NGC 7538 IRS 1 




1.20±0.10 




220+.;$ 


c 70 + O.I5 


5.75-5.83 


5.30-7.00 



, r JV(ortho-HoCO)l 
lQ g[— 277^ L \ 

[cm -2 (km s"i)"i] 
3l2-3l3 4i3-4i4 



13.73] 



14.09±0.05 
14 6 4+0-09 

5+0.03 
-0.01 

is.«±S:i§ 

13.47±g : g? 

14 ^7+0.03 
-0.01 

15.71±0.01 
14.93±0.01 

14.61±0.01 

14 lQ+0.03 
ri.iu_ 04 

14.23±0.01 

13.59±0.05 

14.13±0.01 

14 07+0- 01 

14.03±0.01 
13 50+0.03 

,3 o 7 +0.02 
— 0.01 

14.17±0.01 
13.76±0.01 

-, o 97 +0.07 
1J - ZI -0.09 
10 7 tr+0.03 
10 - '°-0.05 



j+0.05 
-0.06 



^+0.08 
-0.07 



14.18: 

14.64±0.D7 
13.89j 
13.83+° ; °? 
13.64+^ 

14.42±H1 
14.93±0.03 
15.19±0.02 
14 61+0- 04 
14 iq+0.06 
14 ig+0-03 

13 ol+ 014 

14 15+0-05 

14 03+0-05 
1 *- Ul) -0 .09 
14 05+ - 08 
,0 t-7+0.15 

,0 07+0.04 

10,t " -0.03 
14 16+0- 05 

14 04+ - 04 

13 o0+°- 15 
!O.dU_ 24 

13.73+°-?? 



Culled from the literature; in absence of error estimates, a range of 
±50% was tested. Sec i|5.4l for references. 

k Spatial density range found in selected literature. Sec ^5.41 for ref- 
erences. 

c Full LVG modeling was unsuccessful; LTE approximated column 

density quoted. Sec i|5.2l for details on limit derivation. 

^ Total integrated emission used; sec EI5.4.4I ^5.4.91 for details. 



5.4. Comparison to Previous Measurements 

What follows are very brief descriptions of past stud- 
ies of the objects in our sample that were most relevant 
to our analysis. This section also includes accounts of 
any peculiarities in our spectra that resulted in steps 
additional to those described above, notably NGC 1333 
IRS 4A, Orion-KL, IRAS 05338-0624, NGC 2071 IR, and 
DR21(OH). Much has been published on these objects; 
readers seeking additional information should follow the 
citations below. 

5.4.1. W3 IRS 4 

Embedded infrare d source with associated c ompact 
HII Region W3(C) (|Wvnn- Williams et al.lll97l . Our 
spectra found two velocity components at -47.3 and - 
42.5 km s" 1 in the 3i2-3i3 transition, consistent with the 



lm-li i and 2n-2i2 H2CO observations of iDickel et al.1 
(1996), but failed to detect the 4i3-4i4 transition. LVG 
modeling of several H2CO transitions from 211- 365 GHz 
was employed by Mangum fe Woottenl (|1993l ) to find 
T K = 75±f 5 K and log[n(H 2 )] = 5.40+.^ cm" 3 , while 
IHelmich et all (fl99l found T K = 55+.^ K and n(H 2 ) 
= 10 6 cm" 3 in a study using several molecules, including 
four H2CO transitions at around 364 GHz. 

5.4.2. W3(OH) 

Compact HII regio n with a shell structure 
(jDreher fc Welch! Il981h surrounding an embedded 
infrared source with c haracteristics of a young O star 
(iCampbell et al.l 119891) and extensive maser activity 
(e.g. iMoscadelli et all 120101 : Figure [3] of this paper). 
iMauersberger et alj ( 1988l T~detected a warm molecular 
core and red-shifted emission feature using NH3 with 
central velocity ~ -47 km s" 1 . Absorption features, red- 
shifted with respect to the main emission component, 
have also been observed in HC O+ and OH ([Wink et alJ 
Il994t iBaudrv fc Mentenlll995f) . The dominant emission 
and absorption components are detected by both our 
4i3-4i4 and 3i2-3i3 spectra. The J = 3 absorption 
profile exhibits an additional velocity component at 



log[tt(H 2 )] (cm- 3 ) 



Fig. 4. — Comparison of model results for representative examples of the 15 velocity components with transition ratios (Ri) ~ 1-2 that 
were well-fit (left panel) and the remaining 6 components with ratios < 1 whose observed profiles were inconsistent with LVG predictions 
(right panel). Color: |LVG Prediction — Observation! for the integrated intensity ratio Ri (3i2-3i3/4i3-4i4) in units of the measurement 
error a from < 1<t (violet) to > 6<r (red). Contours: |LVG Prediction — Observation] for the 4i3-4i4 brightness temperature T m j in units of 
a\ shaded region corresponds to < 1<t. Observed values listed in the upper right of each panel. Model solution indicated by the intersection 
of the violet and shaded regions. 



47.0kms~ 1 , w hich dominates the 2 and 6 cm H2CO 
observations of iDickel fc Gossl (|1987l ). The absence of 
this component from our J = 4 s pectrum is indicative 
of a f oreground absorbing layer. iMangum fc Woottenl 
(1993) used LVG modeling of several H 2 CO transi- 
tions from 211-365 GHz to find T K = 110±f 5 K and 
log[re( H 2 )] = 5.95^0 01 cm" 3 for the emitting material, 
while IDickel fc Gossj (|1987D suggest that the H 2 CO 
absorption likely arises from material of density ~ 
5xl0 4 cm" 3 assuming a tem p eratu re of 60 K based on 
the work of iGuilloteau et all (|1983f >. They note that 60 
K likely represents an upper limit to the temperature 
of the H 2 CO absorption region and also test 40 K 
to demonstrate the temperature dependency of their 
procedure. 

5.4.3. L1448 IRS 3B 

3A and 3B are Class sources separated by 6" that 
share a common envelope in a potential protobinary sys- 
tem (|Barsonv et al.|[T998f) . of which 3B dominates at mm 
wavelengths dTerebev fc Padgettl 119971 ) H 1 3 CO+ and 
N 2 H + observations bv lVolgenau et al.l (|2006l ) indicate a 
systemic velocity between 3.4-5.8 km s" 1 , consistent with 
our 3i2-3i3 observation (Ulsr = 4.8 km s" 1 ) . The 4i 3 - 
4i4 transition was not detected. iMaret et al.l (|2004[ ) em- 
ployed LVG modeling of several H 2 CO transitions from 
141-364 GHz to find T K = 90 K and n(H 2 ) = 10 5 cm" 3 . 

5.4.4. NGC 1333 IRAS 4A and B 

Premier examples of embedded, low mass star forma- 
tion. Separated by ~ 30", 4A and B a re the bright- 
est o f three sources in the IRAS 4 core (jSandell et al.l 
Il991f ) and each have independently been resolved into 
binary systems (|Loonevi 119981 : lLoonev et al.l l2000h . P 
Cygni profiles, indicative of infall and characterized by 
a blueshifted emission component and redshifted absorp- 
tion, have been detected toward both cores in interfero- 



metric observations of H 2 CO (3i 2 -2n, 226 GHz) and CS 
(|Di Francesco et aill2001l) . This profile is eschewed in all 
but our 3i 2 -3i3 observations of IRAS 4A, which indicates 
an absorption component at 7.5 km s" 1 with two emission 
components at ~ 7.0 km s" 1 being detected in both tran- 
sitions. The broader, wing component of the emission 
can be attributed to outflow material. The absorption 
feature in the 3i 2 -3i3 spectrum of IRAS 4A has the ef- 
fect of bisecting the Gaussian profiles of both emission 
components, eliminating the high velocity component of 
the profile detected in the 4i3-4i4 transition but with- 
out bringing the net J = 3 profile below the baseline. 
Because of the blending in this spectrum, the total inte- 
grated emission over the FWZI of the entire line profile 
for IRAS 4A has been used for our analysis with the inte- 
grated intensity from the Gaussian fit to the absorption 
component added to the total J — 3 intensity. 

The observations toward IRAS 4B were well fit by sin- 
gle Gaussians between 6.9 and 7.1 km s" 1 , in agreement 
with the average syste mic velocity of 7.0 km s" 1 fo und for 
the IRAS 4 group by iDi Francesco et all (|2001fh How- 
ever, the 3i 2 -3i3 emission was anomalously low, result- 
ing in a transition ratio that could not be fit by the 
LVG model. Based on the aforementioned detection of 
H 2 CO absorption toward 4B within the velocity range 
over which our emission was observed, it is possible that 
the J = 3 emission is being subtracted by an absorption 
component not visually reflected in our line profile. This 
possibility is discussed in detail in §6.31 

Using LVG modeling of several H 2 CO transitions from 
141-364 GHz, IMaret et all ()2004P ) derived the following 
properties: 4A-T K = 50 K, n(H 2 ) = 3xl0 5 cm" 3 ; 4B -T K 
= 80 K, n(H 2 ) = 3xl0 5 cm" 3 . Bla ke et alj (|l995l) ana- 
lyzed a combination of H 2 CO (several transitions from 
218-365 GHz) and CS line ratios to find T K = 20-40 K 
and n(H 2 ) = 2xl0 6 cm" 3 for both core components, and 
T K = 70-100 K and n(H 2 ) = 5xl0 6 cm" 3 for the wing 
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material near IRAS 4A. 



5.4.5. L1551 IRS 5 

Low mass star forming site suggested to be a binary 
system separated by ~ 0.3" (jCohen et al.1 fl984) . Sin- 
gle Gaussian emission profiles were detected in both 
transitions with a central velocity typical of the region 
at ~ 6.4 km s" 1 . Mo riartv-Schieven et al.l (|1995f ) used 
H 2 CO (3o3-2 02 and 322-221, 218 GHz) and CS observa- 
tions to derive a kinetic temperat ure of ~ 40 K and a 
spat ial density of 6.9 x lO 6 c m" 3 . iRoberts et all 1)20021 ) 
and lRoberts fc Mil lar (2007) examined transition ratios 
of several species, including the 2n-li (150 GHz) and 
5i4-5i5 (72 GHz) transitions of H2CO, and quote a bolo- 
metric temperature of 97 K for the region. Because the 

TABLE 5 

LVG Results: Velocity Components Detected in One Transition 



T K a tog [ jWaM ]b log[n(H 2 )] 



Source (K) [cm 2 (km s x ) 1 ] (cm" 3 ) 



Emission Components 



L1448 IRS 3B 90 13.17 >5.32 

OMC-2 IRS 4^ ing 85 14.27 4.93<n<6.66 

Cop A East wing 80 13.23 >5.28 

NGC 7538 IRS 9 60 13.39 >5.52 



Absorption Components 

W3 IRS 4 poak 75 13.65 <5.50 

W3 IRS 4 wing 75 13.17 <5.36 

W3(OH) 60 14.22 <5.13 

NGC 1333 IRAS 4A 50 13.71 <5.52 

NGC 2024 FIR 5 25 13.99 <5.58 

NGC 2024 FIR 6 25 13.21 <5.84 

NGC 2071 IR 80 13.48 <5.38 

NGC 7538 IRS 1 25 14.43 <5.17 



H2CO transitions employed by iMoriartv- Schicvc n et al.l 
(1995) are biased toward temperatures < 50 K and the 
LVG prediction of our observed transition ratio sup- 
ports a higher kinetic temperature, we have assumed 
100 K for our analysis while also testing temperatures 
of 40 and 150 K. Additional H2CO measurements have 
been conduc t ed fo r L1551 for the I10-I11 transition by 
lArava et ail (|2006l) . The spatial density of IRS 5 has 
been estimated using several methods. iButner et al.l 
(fl99l found an average volume density of 9x10 cm" 3 
within an angul ar radius of 13'/ 8 by modeling far-infrared 
emission, while iFuller et all (|1995l) estimate 6xl0 5 c m" 3 
using C 17 emission. iMoriartv-Schieven et al.l (|1995[ ) de- 
rived the significantly higher density of llxlO 6 cm" 3 us- 
ing transition ratios of the CS molecule. 



a Culled from the literature. See EI5.4I for references, 
k LTE approximation. 

c For this source only, the 3er detection limit could be used to deter- 
mine upper and lower limits. 

5.4.6. Onon-KL 

The best-studied re gion of massive star for mation to 
date (see review by iGenzel k, Stutzkil I1989D and ex- 
tremely bright in a plethora of chemical species (bright- 
ness temperatures from Orion-KL are > 5x that of any 
other source in our study). This is the only source in 
our sample for which previous observ ations of the 3i2-3i3 
transition of H 2 C O hav e been made (IWilson et al.lll980t 
iMvers fc Buxton] Il980l iBastien et al.lll985D . Outflows, 
shocks, and turbulence arising from newly formed stars 
in the region have led to a complex veloci ty structure 
comp rised of several distinct components (|Blake et al.l 
1987), at least three of which are captured in our ob- 
servations: the hot core ( Vlsr ~ 6 km s , Ay > 10 km 
s" 1 , Tk ~ 300 K), compact ridge ( Vlsr ~ 8 km s" 1 , Av 
~ 4 km s" 1 , Tk ~ 135 K), and extended ridge ( Vlsr ~ 
9-10 km s" 1 , Ay ~ 4 km s" 1 , T K ~ 135 K), each with 
spatial density n(R 2 ) > 10 5 cm" 3 (|Mangum et al.l 119931: 
iMangum fc Wootterj 119931 ) . 

Our J = 4 spectrum exhibits an anomalously intense 



velocity component at Vlsr ~ 10.3 km s , which also 
appears as a much smaller contribution in the J = 3 
transition. A feature a t this velocity exists 10" north of 
our observed position ((Mangum et alJ[T990( ). leading us 
to suspect a pointing error. Pointing was checked prior 
to observing, and the Orion-KL scan was followed by 
observations of OMC-2 IRS 4 and NGC 2024. If a point- 
ing error is to be blamed, it should manifest itself in the 
subsequent observations of OMC-2 IRS 4 and NGC 2024, 
but we could find no evidence for this. The possibilities 
of an unidentified line or rest frequency error were also 
ruled out. This spectrum is the average of only two (mu- 
tually consistent) scans conducted during a single run. 
We suggest that this emission arises from the extended 
ridge, which is typically observed around Vlsr ~ 9 km 
s" 1 , but occasionally as high as 10 km s" 1 . This inter- 
pretation, and the physical parameters derived from it, 
should be applied cautiously for the reasons explained 
above. Because of the disparity between the relative in- 
tensity of the 10.3 km s" 1 feature in each transition, only 
the transition ratio for the hot core component could be 
fit by the LVG model, and the extraordinarily high in- 
tensities precluded the determination of useful density 
limits for the other components through the procedure 
described in N5.21 
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5.4.7. Orion-S 

Broad SiO emission coupled with relatively weak SO2 
and CH3OH emission indicates the presence of an en- 
ergetic outflow in its earliest phase, suggesting that 
Orion-S is one the you ngest stellar objects in the region 
(|McMullin et al.lll993l) . Emission from ^ LSR - 5.5-7.5 
km s" 1 has been obser ved in numerous den se gas trac- 
ers, inclu ding H9CO bvlBastien et al.l (|1985D (2 n -2 12 , 14 
GHz) and Ma ngum et al.l (]1990f ) (several transitions from 
218-291 GHz). Additionally, H 2 CO abso rption has been 
detect ed for the I10-I11 transition by Uohnston et al.l 
(119831) over the velo city range 3.7-9.8 km s" 1 , which 
Mangum et al.l ([19931 ) also detect and attribute to a lower 
density region just north of the emission features. Us- 
ing NH 3 data and the 2n-2i 2 o bservations of H 2 CO 
that would later b e presented by iBastien et al.l |l985), 
iBatrla et al.l (|1983l) found a kinetic temperature of 100 
K for the 6.5 km s" 1 component and a density of 4xl0 6 
cm" 3 , significantly lower than our result. However, they 
also detect a separate velocity component within our de- 
tected velocity range at 7.4 km s" 1 for which they derive 
a density of 6xl0 7 cm" 3 , significantly higher than our 
result. 

We observe both the J = 4 and J = 3 transitions 
in emission with peaks at ~ 6.4 km s . Both spectra 
also exhibit broad wing components likely arising from 
the outflow material. The transition ratio for the wing 
emission was anomalously low and thus the density could 
not be fully constrained by the LVG model. Given that 
H 2 CO absorption has been observed over this velocity 
range, it is possible that the 3i 2 -3i3 emission is being 
partially absorbed without an absorption component be- 
ing visually reflected in the line profile. This possibility 
is discussed in detail in £16.31 



TABLE 6 

LVG Results: Single Transition Observations 



Source 


T K a 
(K) 


log rJV(ortho-H 2 CO)j b 

[cm" 2 (km s" 1 )" 1 ] 
LVG LTE 


Orion-N 
NGC 2024 
G34.26+0.15 
S68N 
W51M 


35 
95 
160 
35 
150 


14 91+0 04 
i4 - zi -0.02 

13.81+°;? 

13.80t°i! 
13.88tg;» 
14.03+°;" 


13.95±0.02 
13.95t°;i° 
14.201°;^ 

13.58lg:i2 
14.401°,;!? 



a Culled from the literature. Sec i|5.4| for references. 
' LVG column density derived assuming n(H2) — 10 61j cm" 3 , the av- 
erage density from Table \4\ LTE results from Table [2] provided for 
comparison. 

5.4.8. OMC-2 IRS 4 

Infrared source composed of two objec ts designated 4N 
and 4 S, which are separated by ~ 4" (jPendleton et al.l 
1986) a nd associated with th e extended (29"xl3") source 
FIR 3 (jMezger et al.l fl990). Our observations exhibit 
Gaussian profiles consistent with the literature and 
peaked at ~ 11.5 km s" 1 . The J = 3 profile includes 
an additional wing component that can be attributed 
to an outflow in the region. Previous LVG modeling of 



several H 2 CO transitions from 211-36 5 GHz was con- 
ducted bv iMangum fc WootterJ (|1993l ). who found T K 
= 85 K an d log[n(H 2 )] = 5.65 cm" 3 , nearly identical to 
our result. iMezger et al.l (|1990f ) estimated the density of 
OMC-2 over a 50"x50" region to be 1.9xl0 6 cm" 3 using 
inferred mass and source size arguments. 

5.4.9. IRAS 05338-0624 

Young stellar object a ssociated with L16 4 1-N, a clus- 
ter of infrared sources ([Strom et al.l 119891: IChen et al.l 
Il993[ ) found to be mainly low - mass, pre-main seq uence 
stars (|Hodapp fc Deand 119931 ). IChen et al.l (|1996l) used 
an analysis of CS observations to estimate the den- 
sity of the mo l ecular core to be ~ 10 6 cm" 3 , while 
iMcMullin et al.l ([19941 ) also used CS to find a density 
of 2xl0 5 cm" 3 assuming a temperature of 42 K de- 
rived from dust continuum obse rvations (| Walker et al.l 
119900 - iStanke fc Williams! (|2007| ) uncovered higher tem- 
peratures (~ 150 K) within a 1'.'4 region using CH 3 CN 
observations. Since the dust continuum represents, at 
best, a lower limit to the kinetic te mperature of the gas 
and th e high temperatures found by Stankc fc Williams! 
(|200l are constrained to an area much smaller than our 
beam size, we have used an average of 95 K for our anal- 
ysis and tested both extremes. 

Both of our H 2 CO spectra exhibit a two-component 
structure with a central peak at Flsr ~ 7.1 km s" 1 and 
a shoulder profile centered around Flsr ~ 9.6 km s . 
The shoulder profile is attributed to the red lobe of the 
outflow material and is consistent with the spec tra of sev- 
eral species observed by [McMunhl£j£aI] (J1994]) ■ Because 
the shoulder profile is very weakly detected in the J = 4 
profile, our confidence in the Gaussian fitting routine's 
ability to reliably separate the two components was low, 
so the total integrated emission over the FWZI of the 
entire line profile was used in our analysis. 

5.4.10. NGC 2024 FIR 5 and 6 

Star forming region containing a string of dense cores, 
FIR 1-7, embedde d in an extended molecular ridge 
(|Mezger et al.l 1992). Both transitions were detected to- 
ward FIR 5 and 6, and an additional 4i3-4i4 measure- 
ment was made for the ridge material (approximately 
midway between FIR 4 and 5). The following prop- 
erties have been derived from LVG modeling of sev- 
eral H 2 CO transitions (211-365, 632 GHz): FIR 5- 
T K = 160±f K, n(H 2 ) = l±0.5xl0 6 cm" 3 ; FIR 6- 
T K = 40+iif K, n(H 2 ) = 2±0.5xl0 6 cm" 3 ; Ridge- T K 
= 95± 3 2 ° K ([Mangum fc Woottenl [19931 : IMangum etHI 
[19991: IWatanabe fc Mitcheill[2008h . 

The absorption feature exhibited by both of our 3i 2 - 
3i3 observations is attributed to a cool (20-30 K) fore- 
ground layer, which has been observed in the lio-ln and 
2n-2i 2 transitions of H 2 CO and found to have n(H 2 ) 
= 10 4 - 9 cm" 3 (IHenkel et al.l[l98l ICrvrFcher et al.1 [19861) . 
The physical parameters of NGC 2024 have been ex- 
plored using a variety of other tra cers as well, recently 
CO bv lEmprechtinger et al.l ([20091 ) . who found the bulk 
of the material to be characterized by Tk ~ 75 K and 
n(H 2 ) ~ 9xl0 5 cm" 3 . 

5.4.11. NGC 2071 IR 
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Cluster of infra red sources spanning ~ 30" 
(|Persson et al.l 119811 ); our observations are centered 
on IRS 1, but IRS 2 and 3 are also being sampled. 
Both spectra indicate a central velocity consistent with 
the literature at ~ 9.7 km s" 1 . An additional shoulder 
component arising from dense gas in the outflow was 
also detected in both transitions at ~ 8.3 km s" 1 , and 
an absorption component was observed in the J = 3 
profile at ~ 16.2 km s" 1 . The 16 km s" 1 feature has been 
pr eviously observed in emis sion for the CS J = 1- line 
bv lTakano et al.l (| 19841 ) and lKitamura et al.l (|1992t ). 

Previous H2CO mea surements of several t ransit ions 
from 211-365 GHz by IMangum fc Woottenl (fl993h es- 
timate T K = 80 K and n(H 2 ) = 10 6 cm" 3 . iTauber et all 
(1988) used the :>> v . 2 I2 (211 GHz) and 3i 2 -2n (225 GHz) 
transitions of H 2 CO to estimate the density in the region 
and found n(H 2 ) ~ 3x 10 5 cm" 3 . A microturbulent model 
of CS and C 34 S emission from lZhou etHI (fl990l ) yielded 
a best fit to the density at 3xl0 6 cm" 3 with the density 
of the outflow emission being > 10 5 cm" 3 . Our LVG anal- 
ysis yielded a significantly higher density for the outflow 
component than for the central emission peak. It is pos- 
sible that our two-component Gaussian fit to the spectra 
is improperly separating the components. An analysis of 
the total integrated emission over the FWZI of the entire 
line profile yields ?i(H 2 ) = 10 6 00 cm" 3 . 

5.4.12. S255N 

Massive star-forming region also known as S255 FIR 1 
or G192.60-MM1 that lies at one end of an extende d 
molecular ridge opposite S255IR (Hcv er et al.l I1989I) . 
Three compact cores, SMA 1-3, were resolved in th e 
1.3 mm continuum maps of iCyganowski et"aTl (|2007[ ). 
who also combined an analysis of H 2 CO transition ra- 
tios (3 03 -2o2, 3 22 -2 2 i, and 3 2 i-2 20 ; 218-219 GHz) and a 
spectral energy distribution (SED) model to estimate a 
temperature range of 40-100 K and densities between 
3-16 xlO 6 cm" 3 . The H 2 CO observations presented in 
ICyganowski et all ()2007l ) are centered on either side of 
SMAl (NE and SW) and indicate velocity components 
at V^lsr ~ 6-9 km s" 1 and 12.1 km s" 1 , while our spectra, 
centered between the two, show a single component at ~ 
8.9 km s" 1 . 

5.4.13. G34.26+0.15 

Ultra compact (UC) HII region with an associated hot 
core that has become a prototypical example of cometary 
morp hology (jWood fe Churchwellll 19891 : Ivan Buren et al.l 
119901 ). The hot molecular gas (80-175 K) is suggested to 
be the outer layer of a massive, cool core that is be- 
ing externally heate d by the UC HII region, from which 
it is offset by - 2" (|Heaton et al.lll989|- IWatt fe Mundvl 
1999). Single di s h CH 3CN observations conducted by 
Churchwell et al.l (|1992D suggest gas temperatures of 166 
K an d densities > 10 5 cm" 3 , while iMookeriea et al.l 
(2007) also found a temperature of 160 K using the 
brightness temperatures of several optically thick lines. 

5.4.14. S68N 

Deeply embedded source within the Serpens molecular 
cloud containing a proto star with an as sociated outflow 
(Wolf -Chase et alj I1998D . suggested by iMcMullin et al.l 
2000) to be an example of a "cool core," with 



properties intermedia te to cold/warm condensations. 
IMcMullin et al.l (|2000|) used relative intensities of four 
H 2 CO transitions [l O i-0 O0 (73 GHz), 2 02 -l i (146 GHz), 
3o3-2o2 (218 GHz), and 3 22 -2i 8 (218 GHz)] to derive a 
kinetic temperature range of 35-70 K but report flat- 
topped line profiles suggesting an overestimation due to 
optical depth effects and adopt a n estim ate of 35 K for 
their measurements. iHurt et ail (|1996l ) also examined 
four H 2 CO transitions [3 03 -2 02 (218 GHz), 3 22 -2i 8 (218 
GHz), 5o5-4 04 (363 GHz), 5 23 -4 22 (365 GHz)] to deter- 
mine a higher temperature of around 75 K. Both studies 
also re port spatial density estimates withlMcMullin ^et al.l 
(|2000l) finding n(H 2 ) =0.4-1.2xl0 6 cm" 3 using H 13 CO+, 
SiO, and DCN, while IHurt et~aT1 (fl996l ) found 2.5xl0 6 
cm" 3 using H 2 CO. 

5.4.15. W51M 

Do minant r e gion o f the massive star- forming site W51, 
which Martin (|1972| ) revealed to consist of eight distinct 
components at centimeter wavelengths. Our observa- 
tions are centere d on W51e, itself di vided into four UC 
HII regions el-4 (jGaume et al.i ri993). whose association 
with dense m olecular cores and maser activity is known 
as W51-Main (jZhang fe Hd[T997l ). iRemijan et all (pOOl 
used CH3CN to find temperatures of 123 and 153 toward 
el and e2, re spectively, and a de nsity of 5xl0 5 cm" 3 for 
both sources. iZhang fe Hoi |l997) derive higher densities 
of 2-3 xlO 6 cm" 3 for the region using NH3. 

5.4.16. DR21(OH) 

Star-forming region in an early phase of its evolution 
which has yet to see substantial ionization of the molecu- 
lar material surround ing newly-formed, massive B stars. 
IMangum et al.l (|1992t ) detected four principal condensa- 
tions labelled M, N, W, and S. Our observations are 
centered on DR21(OH)-M, which is subdivided into two 
components separated by ~ 8", MM1 and MM2. De- 
spite their proximity, MM1 and MM2 differ substantially. 
MM1 is hot (~ 160 K), dense (~ 10 6 cm" 3 ), and moder- 
ately luminous (~ 10 4 Lq), while MM2 is cooler (~ 30 
K), denser (~ 1 7 cm " 3 ), and less luminous (~ 10 3 Lg ) 
(IMangum et al.lll991l [l992t IMangum fc Wootten1ll993f) . 
Our J — 3 spectrum indicates distinct peaks at -4.5 
and -1.3 km s" 1 , attributed to MM1 and MM2, respec- 
tively. The J = 4 spectrum is severely blended and de- 
mands a surprisingly high contribution from the -1.3 km 
s" 1 (MM2) component. Since this region has not been 
mapped at these frequencies and the spectra were so 
highly blended, due caution should be exercised when 
assessing our results. 

5.4.17. Cep A East 

Massive star forming region known to consist of 16 
compact (~ 1") components clustered within a 25" 
radius and aligned in an inverted Y-shaped structure 
(|Garav et al.l fl996V Our observations are centered on 
the continuum source HW3, ~ 3" south of the domi- 
nant source HW2. HW2 is well within our beam and 
is associated with a pr omising candidate for the detec- 
tion of a massive disk (| Torrefies et al.|[l996l : iPatel et ail 
2005), though other studies suggest that this elongated 
molecular structure is explained by the superposition of 
at least three hot cores (|Comito et al.1 2007; Broga rTet al.1 
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2007) . Choosing an appropriate kinetic temperature esti- 
mate was difficult for this region because of the inclusion 
of several compact sources within our beam. We have 
as sumed a rang e of 40-100 K based on the estimates 
of lBrogan et~aTT p007t ) and lCodellaeTall pOOll) for the 
ambie nt cloud velocity of ~ 10.5 km s" 1 . Codell a~et al.l 
(2005) also find densities between 2-6 xlO 6 cm" 3 in the 
region using SO and SiO. 

5.4.18. NGC 7538 IRS 1 

The brighte st of three compact sources discovered by 
iMartinl (|1973f ). which contains a zero-a ge main sequence 
star of spectral type earlier than 07.5 (jLugo et al.ll2004D 
embedded in a molecular cl oud with an inner sh ell struc- 
ture surrounding the star (jPratap et al.l 119971 ). A sys- 
temic velocity of around -57 km s" 1 was observed in 
several transitions of H 2 CO between 365-470 GHz by 
Ivan der Tak et all (|2000l ) with an absorptio n feature near 
-60 k m s" 1 having been detected in NH 3 (|Wilson et al.l 
fl983h . both of which are found in our spectra. The 
region is also notable f or the rare occurrence of 6 cm 
(lio- lii) H2CO masers (jRots et al.lll98lHHoffman et al.l 
l2003f ). The absorption componen t we observe was also 
detected by iHoffman etaLl (|2003ft in the 2cm (2ii-2 12 ) 
transition of H2CO. 

A temperature of ~ 220 K was found using NH3 by 
IMauersberger et aL (1988), while a cooler temperature 
of 176 K was found bv lMitchell et~aTl ffl990h using 13 CO 
and a warmer one of 245 K bv IQiu et al.l (I2TTI using 
CH3CN. The Mitchell study also indicates a cold gas 
component of 25 K from which the absorption compo- 
nent in our J = 3 spec t rum likely arises. From 13 CO 
observations, IQiu et all (|2011l ) estimate the density to 
be ~ 10 7 c m" 3 for a region abo ut half the size of our 
beam, while [Mitchell et al.l ()1990[ ) indicate that the den- 
sity is > 10 6 cm" 3 . IHoffman et al.l (I2003D note that to 



excite the lio-ln H 2 CO masers in the region, lower den- 
sities in the range of 6-16 xlO 4 cm" 3 are required. The 
masers arise from a region directly in front IRS 1 not 
coincident with the hot core region examined by the pre- 
viously mentioned authors. 

Since our observations are of higher excitation transi- 
tions than the H2CO masers in the Hoffman study, it is 
likely that we are sampling material from both the hot 
core and maser regions, and thus the intermediate den- 
sity we find is justified. However, this notion also calls 
into question the validity of our kinetic temperature as- 
sumption, so it should be noted that temperatures < 120 
K would result in densities > 10 6 in the LVG approxi- 
mation. That said, due to the excitation requirements of 
the 3i2-3i3 and 4i3-4i4 transitio ns compared to that o f 
the J = 1 JT-doublet observed bv IHoffman et al.l (|2003j) . 
it is more likely that our measurements are biased toward 
the higher temperatures of the hot core. Fortunately, the 
H2CO LVG models are relatively independent of Tk for 
temperatures > 100 K (see 



5.4.19. NGC 7538 IRS 9 

Deeply embedded cold IR source ~ 50" south of 
IRS 1 that is associate d with a large reflection neb- 
ula (|Werner et al.l I1979D . The radial velocity of the 
IRS 9 cloud core from H 13 CN CS, and H 2 CO data is 
about -57 km s" 1 dSandell et alj|200a : Ivan der Tak et all 



2000), consistent with our observations. Single-dish 
JCMT H 2 CO (3 22 -2 2 i and 3 2 i-2 20 , 2 1 8 GHz) spectra sug - 
gest a gas temperature of > 60 K (|Sandell et al.ll200"5T) . 
while HCO+, H 13 CO + , and 1 3 CO observations suggest a 
lower temperature o f ~ 30 K ( Hascgawa fc Mitchelil995l: 
iMitchell et all 119901 ) . IMitchell et all (|1990l) alsodetect 
the presence of warm gas (180 K) and suggest that the 
density of this material is > 10 6 cm" 3 . 

6. LIMITATIONS OF H 2 CO J = 3/J = 4 if -DOUBLET 
DENSITOMETRY 

Our study is concerned primarily with the appli- 
cation of a previously unused densitometry technique 
and, as such, well-studied objects were chosen to as- 
sess the method's performance. This discussion is fo- 
cused on limitations to the technique that require elab- 
oration. These effects may also help explain the prob- 
lem of the anomalously low transition ratios [J T m b(3\2 — 
3i3)df/ / r m j,(4i3— 4u)di'] that precluded LVG modeling 
in a few sources (see Table UJ. 

TABLE 7 

Temperature Dependence of LVG Density Derivation 



T k 


log[n(H 2 )f 


% Difference 


(K) 


(cm -3 ) 


Over Tj( Range 


30-50 


7.72-6.71 


923 


50-100 


6.71-6.12 


290 


100-150 


6.12-5.91 


59.2 


150-200 


5.91-5.81 


26.2 


200-250 


5.81-5.75 


15.0 


250-300 


5.75-5.69 


15.0 



Calculated for tost case NGC 7538 IRS 1. 



6.1. LVG Model Dependence on Kinetic Temperature 

The LVG model for the 3i2-3i3 and 4i3-4i4 transitions 
is somewhat dependent on kinetic temperature (Tk)- As 
previously mentioned in H5.ll this meant that a careful 
selection of appropriate kinetic temperatures from previ- 
ous studies was necessary. Wherever possible, tempera- 
tures derived from H2CO measurements were used to en- 
sure coupling to the gas traced by our observations. Esti- 
mates taken from analyses of other dense gas tracers such 
as NH 3 or CH 3 CN are otherwise preferable. That said, 
any time a temperature derived from measurements of 
disparate molecules (or even H2CO measurements of sig- 
nificantly different excitation requirements) is adopted, 
the question of whether or not this estimate can be asso- 
ciated with the gas sampled by our beam must be raised. 
In at least one case (NGC 7538 IRS 1, 35.4.18P there ap- 
pears to be a distinct possibility that the temperature we 
adopted may not be wholly appropriate for the material 
traced by our observations. 

In Table |H we elected to include the spatial densities 
derived for a range of kinetic temperatures based on the 
error estimates in our assumed values. This is because 
the temperature dependence of our H2CO LVG models 
is not linear and itself depends on the range of physical 
parameters being studied. Table [7] provides a summa- 
tion of this dependence by dividing the range of kinetic 
temperatures found in our sample into 6 groups. The 
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test case of NGC 7538 IRS 1 was chosen for its mod- 
erate transition ratio (Ri) of 1.20. From the table it is 
clear that the dependence on Tk decreases with an in- 
crease in the assumed kinetic temperature. In the most 
volatile range (30-50 K) the measured density can fluc- 
tuate by upwards of an order of magnitude between the 
extremes. Fortunately, very few of our sources fall in this 
range. At high temperatures (> 200 K) the dependence 
on Tk is fairly constant. In this regime, uncertainties in 
the derived density due to errors in the adopted kinetic 
temperature are generally below the uncertainty imposed 
by the collisional excitation rates between H2CO and H2 
(20%; see gET} . 

Figure [5] provides a graphical representation of the 
H2CO LVG model's dependence on Tk with plots identi- 
cal to that of Figure|4]for a source of high and low kinetic 
temperature. Note how at high temperatures (upper 
panel), the transition ratio and brightness temperature 
contours approach orthogonality. In effect, this means 
that as the kinetic temperature varies and the locations 
of the curves change with respect to each other, the inter- 
section point between them is moved by smaller amounts 
at higher temperatures. For Tk > 200 K, the brightness 
temperature (solid-line) contours become roughly hori- 
zontal and the solution point is therefore only affected 
by the lateral motion of the transition ratio contours. 
This is why a variation in Tk from 200-250 K produces 
the same variation in spatial density as 250-300 K (see 
Table El). 

It should also be noted that, depending on the observed 
J = 3/J = 4 integrated intensity (line) ratio, the kinetic 
temperature assumption can be the deciding factor in 
whether or not the LVG model can successfully match 
an observation. At high temperatures, the higher excita- 
tion J = 4 transition can start to enjoy a greater relative 
population than the J = 3, and Rf < 1 can reasonably 
be expected. In Table 0] there are 3 sources, L1551 IRS 
5, Orion-S, and DR 21(OH)-MMl, with transition ratios 
on order ~ 0.9. This was not a problem for the assumed 
temperature, but at the lower temperature limit, the den- 
sity is indicated to be > 10 8 (column 5). This is to say 
that the LVG prediction for the density is unbound at 
such low temperatures. In terms of Figure [SJ for kinetic 
temperatures below a given value, the transition ratio 
(color) and brightness temperature (solid-line) contours 
no longer intersect (as in Cep A East from Figured]). 

Given this discussion and the uncertainties typical in 
estimates of kinetic temperature, we suggest that H2CO 
3i2-3i3/4i3-4i4 densitometry is best suited to objects 
with Tk > 100 K. We have demonstrated that this tech- 
nique is also viable at lower kinetic temperatures but 
with the caveat that its results become increasingly sen- 
sitive to uncertainties in the assumed value of Tk- The 
reason for this effect can likely be attributed to the rel- 
ative excitation requirements of the J = 3 and 4 transi- 
tions, which are discussed in the following section. 

6.2. H2CO 3i%-3i3 and 4i3~4i4 Relative Excitation 
Requirements 

In addition to the kinetic temperature of an object, it 
is important to consider the relative excitation require- 
ments for each transition. Being of lower excitation, the 
3x2-3x3 transition can be excited by temperatures and 
densities lower than that required by the 4x3-4x4 transi- 



< g >6 




4.5 5.5 6.5 



log[n(H2)] (or') 

Fig. 5. — Comparison of model results for sources of high (up- 
per panel) and low (lower panel) kinetic temperature. Contours 
and parameters are as in Fig. [4] Note how at high Tk, contours 
approach orthogonality, resulting in less variation in the solution 
point due to changes in Tk (alternatively, less dependence on Tk)- 

tion. This means that the J = 3 transition is predictably 
observed to be more intense than the J = 4 in the general 
case. What complicates the matter is that the two tran- 
sitions are affected differently by the collisional pumping 
mechanism described in |J21 In effect, this means that 
the 4x3-4i4 transition falls into absorption at a slightly 
higher spatial density than the 3x2-3x3- Therefore a small 
band of parameter space, briefly described in §5.2[ exists 
for which the J = 4 transition is in absorption while the 
J = 3 is in emission. 

It may be counterintuitive, then, that only one absorp- 
tion component was detected in the 4x3-4x4 transition 
while 9 were observed in the 3x2-3x3- This is explained 
by two effects. The first is that the system temperatures 
for the J = 4 (Q-band) observations were substantially 
higher than the J = 3 (Ka-band) , resulting in much nois- 
ier 4x3-4x4 spectra. Looking at the spectra in Figure [TJ 
the majority of absorption components detected in the 
3x2-3x3 transition are below the noise of the 4x3-4x4 spec- 
tra. The second, and perhaps more significant, effect is 
that the J = 3 beam is substantially larger than that of 
the J — 4, meaning that the material being sampled is a 
bit different for each transition. This topic is addressed 
in AO 

Figure [6] displays the excitation temperatures of the 
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H 2 C0 K-Doublet Excitation Temperature 
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Fig. 6. — H2CO J = 3 and 4 if -doublet excitation temperatures 
as a function of molecular hydrogen density for an ortho-H2CO 
column density per velocity gradient roughly equivalent to the av- 
erage value observed in our sample. Note that the possibility for 
absorption is much smaller at high kinetic temperatures than it is 
for low. Values predicted via LVG modeling. 

J = 3 and 4 transitions as a function of spatial density 
for an ortho-H 2 CO column density roughly equivalent to 
the average of that observed in our sample. From this, we 
can see the effect of the collisional pumping mechanism 
described in ij2j which cools the excitation temperatures 
of the J < 5 if -doublet transitions to lower than 2.73 K 
for a given set of physical conditions, allowing them to 
absorb radiation from the cosmic microwave background. 
Curves for kinetic temperatures of 50, 150, and 200 K are 
shown. Note how at high kinetic temperatures, the ex- 
citation temperatures for both transitions match more 
closely and only barely fall below the 2.73 K level. This 
means that at high kinetic temperatures, the possibility 
of observing absorption is small and the difference be- 
tween the excitation requirements of each transition is 
of less consequence. This is likely the reason why the 
3i2-3i3/4i3-4i4 densitometry technique is best suited for 
temperatures > 100 K. 

6.3. The Effects of Spatial Extent 

Because the beam sizes differ by (6>fc(29 GHz) — 
b {4&GHz))/0 b (ASGHz) = 63%, it is important to con- 
sider the distribution of dense gas in each object as the 
larger beam size observations (J = 3) may be sam- 
pling extended structure that biases the resulting line 
(J = 3/J = 4 integrated intensity) ratio. In similar stud- 
ies using the I10-I11 and 2n-2i 2 if -doublets, the beam 
sizes of both observations (153" and 51", respectively, for 
the GBT) can generally be considered much larger than 
the dense gas distribution and it is safe to assume that 
all relevant emission is sampled. One advantage of the 
3i2-3i3 and 4i3-4i4 transitions is that their smaller beam 
sizes (26" and 16", respectively, for the GBT) allow for 
higher resolution and heightened sensitivity to the spa- 
tially compact regions in which stars are known to form. 
However, these beam sizes also approach the spatial ex- 
tent of molecular cores and it is therefore prudent to con- 
sider potential differences in the material being sampled 
by each observation. 

Unfortunately, no mapping measurements of the H2CO 
J = 3 and 4 if -doublet transitions exist. We experi- 
mented with assuming spatial extents derived from other 



dense gas tracers in the literature but ultimately decided 
that inconsistencies in the frequency, chosen tracer, exci- 
tation conditions, and spatial resolution precluded truly 
appropriate spatial extent estimates for our purpose. In- 
stead, source sizes equal to the smaller beam size (16") 
were assumed and the J = 3 observations were scaled for 
beam dilution (see We find this agreeable because 
16" is a reasonable average value for estimates of the 
dense gas distribution in most of our sources, but this is 
surely not appropriate for each application. 

It is also possible for an object's dense gas distribu- 
tion to be significantly smaller than 16", in which case 
the J = 4 observation is also subject to beam dilution. 
To examine this effect, consider the case of L1551 IRS 
5, for which high- resolution imaging of CS suggests the 
dense gas to be distributed over ~ 7" (jTakakuwa et al.1 
2004). Applying the beam dilution correction described 
in £J3] to this case has the effect of boosting the inten- 
sity of both transitions, but comparatively more for the 
J = 3 observation, resulting in an increased transition 
ratio [/ T m ;,(3i2 - 3i 3 )dv/ / T m &(4i 3 - Au)dv\. The sub- 
sequent LVG approximation to the density (10 5 95 cm" 3 ) 
is nearly an order of magnitude lower than that of the 16" 
assumption (10 6 87 cm" 3 ). This decrease is mainly due to 
an increased optical depth in the modelled transitions. 

In several sources, the transition ratio was anomalously 
low (see Table |4]), meaning that the J = 3 transition 
was observed to be unexpectedly less intense than the 
J = 4. This is possible at temperatures and densities 
high enough to populate the higher excitation transition 
significantly more than the lower, but these conditions 
are not met by the sources in question. One possible ex- 
planation for this effect is that the J = 3 transition is 
being subjected to absorption that is not reflected in the 
line profile other than to lower its intensity. The LVG 
model assumes a uniform spatial density, which proves 
to be a reasonable approximation but ignores a density 
gradient that likely exists in molecular cores, with the 
outer regions being comparatively less dense than the 
inner. It is possible that the larger J = 3 beam is sam- 
pling the comparatively cool, diffuse material enveloping 
a molecular core which has resulted in some self absorp- 
tion. Indeed, in a few sources (e.g. NGC 1333 IRAS 4B 
and Orion-S, 35.4.41 and 35.4.71 respectively), H 2 CO ab- 
sorption has been detected in the J = 2 or 1 if -doublet 
transitions over the velocity range in question. Mapping 
studies are required to disentangle the varying density 
structure in many of the sources in our sample. 

7. CONCLUSION 

Using observations of the 3i2-3i3 and 4i 3 -4i4 transi- 
tions of H 2 CO, we have successfully constrained the spa- 
tial densities of a sample of galactic star-forming regions. 
Both transitions were observed toward 18 sources with 
relative ease, requiring an average of 17 min of integra- 
tion time and resulting in only 3 nondetections in the 
J = 4 transition. Accurate measurements of the spa- 
tial density [ra(H2)] were made for 13 objects and use- 
ful limits were placed on the remainder using a com- 
bination of Large Velocity Gradient (LVG) and Local 
Thermodynamic Equilibrium (LTE) analyses. Molecular 
hydrogen densities in the range of 10 5 ' 5 — 10 6 5 cm" 3 and 
ortho-formaldehyde column densities per unit line width 
between 10 13 5 and 10 14 5 cm" 2 (km s" 1 )" 1 are found for 
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most sources, in general agreement with previous mea- 
surements. 

Detailed analyses of the advantages and limitations 
to this densitometry technique have also been provided. 
H 2 CO 3i2-3i3/4i3-4i4 densitometry proves to be best 
suited to objects with Tk > 100 K, above which the 
H2CO LVG models become relatively independent of ki- 
netic temperature. Compared with the similarly utilized 
I10-I11 and 2n-2i2 transitions, the J = 3 and 4 K- 
doublets provide higher spatial resolution and sensitivity 
to hot, dense material, which makes them more efficient 
probes of spatial density in molecular cores. However, 
beam widths comparable to the anticipated source sizes 
make source structure considerations important. Since 
mapping measurements have yet to be conducted for 
these transitions, the correlation between the spatial ex- 
tent traced by the J = 3 and 4 if -doublets and that 
of other dense gas tracers is uncertain, making spatial 
assumptions based on past measurements problematic. 



This work serves as a successful proof-of-concept for 
the H 2 CO J = 3/J = 4 if -doublet densitometry tech- 
nique, adding a useful new diagnostic to the study of 
dense molecular environments. 
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